Complex orbital state in manganites 



Ryo Maezono and Naoto Nagaosa 
Department of Applied Physics, University of Tokyo, Bunkyo-ku, Tokyo 113-8656, Japan 

(February 6, 2008) 

Abstract 

The e g -orbital states with complex coefficients of the linear combination 

of x 2 — y 2 and 3z 2 — r 2 are studied for the ferromagnetic state in doped 

manganites. Especially the focus is put on the competition among uniform 

complex, staggered complex, and real orbital states. As the hole-doping x 

increases, the real, the canted complex, and the staggered complex orbital 

states appears successively. Uniform complex state analoguous to Nagaoka 

ferromagnet does not appear. These complex states can be expressed as a 

resonating state among the planer orbitals as the orbital liquid, accompanied 

by no Jahn- Teller distortion. 
71.27.+a, 75.30.-m, 75.30.Et 
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I. INTRODUCTION 



The role of the orbital degrees of freedom has recently attracted considerable interests 
as one of the key to understand the colossal magneto-resistance (CMR) observed in doped 
manganites.§Hll The orbital state of the conduction electrons is described as a linear com- 
bination of two wavef unctions, \x 2 — y 2 ) and \3z 2 — r 2 ), of the degenerate e g orbitals.0 In 
previous studies,0lll the linear combination with only real coefficients (real orbital state) has 
been considered. This is because theories of the orbital ordering have been developed mainly 
to describe the parent compounds of CMR materialsjl!~0 in which the static Jahn- Teller 
deformation is observed.0 Such a deformation stabilizes the real orbital state and it was rea- 
sonable to exclude the linear combination with complex coefficients (complex orbital state). 
Recently the orbital state in doped compounds is studied concerning the properties of CMR 
materialsllHl! Because the static Jahn- Teller distortion disappears in doped compounds JllEiil 
there is no reason to exclude the complex orbital state. Actually such a complex orbital 
state has been recently studied.!!! 2 ! Khomskii0 pointed out that the complex orbital state, 
(\x 2 — y 2 ) ± i \3z 2 — r 2 )) j y/2, provides locally isotropic hopping intensities with the same 
bandwidth as the real orbital state, and might explain the isotropic properties observed in 
CMR compounds. Such a local isotropy cannot be realized with the uniform real orbital 
state.S A staggered ordering is therefore needed to explain the observed isotropic properties 
within the extent of the real orbital ordering^ (Another proposal is the orbital liquid state, 
where the local isotropy is recovered by a quantum resonance between anisotropic orbital 
configurations, \x 2 — y 2 ), \y 2 — z 2 ), and \z 2 — x 2 )§). Based on the analogy to the Nagaoka 
ferromagnetism (F) , Khomskii proposed that the uniform ordering (orbital F) of the com- 
plex orbital state is more stable than the staggered one (orbital AF) with real orbitals.!! 
Takahashi et al. investigated the possible complex orbital ordering, motivated by the anal- 
ogy to the octapole ordering in heavy fermion systems with odd time reversal symmetry.!! 
They found that the staggered ordering of the complex orbital is stable, being contrary to 
Khomskii's uniform one.@ 
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In this paper, we study the competitions among the uniform complex, staggered complex, 
and real orbital states by using a model of CMR compounds taking the strong on-site 
repulsion and the orbital degeneracy into account .0 The complex orbital state is taken as, 
cos(6 l /2) ■ \x 2 — y 2 ) + i sin (0/2) • \3z 2 — r 2 ), and the whole possibility with the continuous 
parameter 9 is examined. With realistic parameters, the complex orbital state is more stable 
than the real one in the moderately doped region (0.25 < x < 0.45). The complex ordering 
changes from the canted one (0.25 < x < 0.35) into the staggered (0.35 < x < 0.45) one 
due to the competition between the orbital superexchange AF and the orbital Nagaoka F. 
The local isotropy is also realized in this complex staggered phase, where the band gap 
due to the doubled period brings about the energy gain exceeding the energy loss due to 
the narrower bandwidth than that of the uniform ordering with isotropic hopping. With 
increasing U/t toward the strong correlation limit, the former gain decreases whereas the 
latter loss increases. The staggered ordering becomes unstable in this limit, where the 
uniform orbital ordering wins. In this case, however, the obtained uniform ordering is not 
the complex one0 but the real one with \x 2 — y 2 ). Though the uniform complex ordering 
becomes more stable than the staggered complex one, it has higher energy than that of the 
real one. In the weak correlation limit, on the other hand, Takahashi et al. found that the 
normal metallic state becomes unstable toward the the staggered complex ordering near the 
quarter filling,!! with increasing U. When the Jahn- Teller coupling is further taken into 
account, however, it is likely that the real orbital state is stabilized, because the energy 
scale of the Jahn- Teller coupling becomes dominating compared with the weak U, prefering 
the real state. 

The Jahn- Teller deformation which couples with the orbital degrees of freedom decreases 
in the complex canted phase and vanishes in the complex staggered phase, being consistent 
with the observed disappearance of the deformation. 00 This complex state can be expressed 
as a resonating state among planer orbitals, as in the orbital liquid picture.^ When the 
resonance occurs with coherent correlations in time and space, the complex orbital ordering 
is obtained, meanwhile that with incoherent one corresponds to the orbital liquid state.@ 
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These can be distinguished by experiments detecting the spatial correlation of the orbital 
symmetry, such as the anomalous X-ray scattering experimentsSH Possibilities of the phase 
separation and broken time-reversal symmetry are also discussed. 



II. RESULTS AND DISCUSSIONS 

We employ the same model as that in the previous report!, 

0-77' {ij) 



+ Js X! St 2 gi'St 2g j + -^on site, (1) 
(ij) 

where 7 [= a(d x 2_ y 2), b(d 3z 2_ r 2)] specifies the orbital and the other notations are standard.0 
The transfer integral t]j depends on the pair of orbitals (7, 7') and the direction of the bond 
(2, The spin operator for the e g electron is defined as S egi — | dl^a^d^p with the 
Pauli matrices a, while the orbital isospin operator is defined as Tj = | Y dl^a^ydiy^B 
Jh is the Hund's coupling between e g and t 2g spins, and Js is the AF coupling between 
nearest neighboring t 2g spins. H on site represents the on-site Coulomb interactions between 
e g electrons. Coulomb interactions induce both the spin and orbital isospin moments, and 
actually H on site can be written as 

H on sitc = -£ (ffi + aS 2 egt ) . (2) 

i 

A parameter set with to = tf b i+i = 0.72 eV, a = 8.1t , and (3 = 6.7to corresponds to the real- 
istic one being relevant to the actual manganese oxidesE In the path-integral quantization, 
we introduce the Stratonovich-Hubbard fields (ps and <pt, representing the spin and orbital 
fluctuations, respectively. With the large values of the electron-electron interactions above, 
both 0s and ipx are almost fully polarized.0 The meanfield theory corresponds to the saddle 
point configuration of ips and 0T- We only consider the possibility of the complex orbital 
state within a F-type spin alignment in the cubic cell. 



We assume the two sublattices, I and 77, with F-, A-, C- and G-type alignment. On 
each site, the orbital is specified as a linear combination of the two degenerate orbital bases, 
\ x<2 ~ V 2 ) an d |3z 2 — r 2 ), as 
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9,<p) = cos- x 2 -y 2 ) + e"^sin- 3z 2 - r 2 ) . (3) 
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(9, if) is the polar angle of the corresponding isospin T. In the limit of the infinite orbital 
polarization, (3 — > oo, the uniform orbital ordering with any \9, if) takes the same bandwidth, 
— (3/2) to- The polar angle (9, if) therefore controls only the dimensionality of the band 
structure to optimize the kinetic energy gain, leaving the bandwidth unchanged. 

Previous studies^!! have focused on the states with T j /e y as the complex orbital states. 
In this paper, we extend the possibility to T lying within yz plane (\9,ip = tt/2), real and pure 
imaginary coefficients) for the complex orbital state, whereas T within zx plane (\9,if = 0)) 
corresponds to the real orbital state. This choice includes \x 2 — y 2 ) and |3^ 2 — r 2 ) as the 
both ends. With finite (3, the generalized orbital canted structure on two sublattices is 
examined. 

Fig. 1 shows the energy values in spin F phase, optimized within the real and the 
complex orbital states, plotted as a function of the hole concentration x (with Js = 0). 
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FIG. 1. Energies of the spin F phase with real and complex orbital states as a function of the 
hole concentration x (with J5 = 0). 



The orbital shape specified by 6 is optimized at each x. The complex orbital state is 
realized in the moderately doped region (0.25 < x < 0.45). The phase diagram as a function 
of x and J$ {AF interaction between t 2g spins) is shown in Fig. 2. 




coated staggered 

FIG. 2. Phase diagram as a function of the hole concentration (x) and the antiferromagnetic 
interaction between t 2g spins (Js). A, C, F, and G specify the spin configuration. 



In the shaded and hatched regions of the spin F phase is realized the complex orbital 
state. The phase boundary depicted with a broken line is that for the real orbital state, 
reported previously!] 
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FIG. 3. The phase diagram assuming the spin F phase with realistic parameters. 

Figure |3| shows the orbital phase diagram assuming the spin F phase as a function of 
x. The orbital ordering changes from the real staggered, the complex canted, the complex 
staggered, and to the real uniform one. Figure [| shows the x-dependence of the orbital 
canting angle, \6n — 6i\, for the real and the complex orbital states. 
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FIG. 4. x-dependence of the orbital canting angle for the real and the complex orbital states 
(with realistic parameters). 

With increasing x, the canting angle once tends to take the orbital F (I — > //), but 
get back to AF again (II — > III). The canted complex state with orbital C is stable at 
x = 0.25 and 0.35. With increasing x, the canted state changes into the staggered one for 
0.35 < x < 0.45 with = —6 U = n/2 (orbital G) as found in ref. 22. We note that this 
staggered state also gives the locally isotropic hopping integrals, t x = \e~ %2 f , t y = |e*T", 
and t z = — |. The uniform complex state0 with 6j = 6jj = n/2 and t x,y,z = — 1/2 has higher 
energy. With further doping (IV with x > 0.5), the orbital F becomes stable again, but 
with real coefficients. 

These results can be understood as follows. The orbital superexchange AF interaction 
Jaf is represented by the shift in the center of mass of the occupied density of states (DOS), 
as represented by the Hamiltonian, 

7 



£fc PeB 
. PeB e k+Q 

with Q = (tt, tt, tt) being the staggered orbital momentum. Therefore Jaf is estimated as 



(4) 



J 



t 2 



t 2 



AF 



(5) 



PeB P(l-X) ' 

which increases as x increases because p e g is the constant (3 times the number of electrons 
(1 — x). The ferromagnetic double exchange interaction Jf for the orbital moments is rep- 
resented by the energy of the doped holes at the top of the occupied band, which depends 
on the bandwidth. The bandwidth is t for the uniform ordering whereas t 2 / '/3 e s for the 
staggered one for t <C P e g and small x. Jp is therefore given as, 

i t 2 \ 



Jv 



t 



X 



'eff. 



(6) 



which represents the relative kinetic energy gain of the orbital F state measuring from that 
of the staggered state. It should be noted here that the notation Jp is rather symbolic, and 
the Hamiltonian is not written as — Jf S« Ti ■ Tj. Based on these considerations, the results 
in Fig. |2] and ||] are interpreted as follows. Here we assume the ferromagnetic spin alignment. 
At x = 0, Jp vanishes meanwhile Jaf is finite, leading to the orbital AF. With small doping, 
Jf becomes finite, leading to the tendency toward the orbital F seen in the region II in Fig. 
[|. (This corresponds to the crossover from the orbital superexchange AF to the orbital 
double exchange (Nagaoka) F with the doping.) To understand the reentrant of the orbital 
AF in the region III, we note that t 2 //3 e s = t 2 / (3 (1 — x) increases as x, which enhances 
Jaf (x) and suppresses Jp. Actually, the difference in the bandwidth of DOS between the 
uniform (orb. F) and the staggered (orb. AF) structures is hardly seen at x = 0.3 in Fig. 
|5], corresponding to Jf ~ 0. 
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FIG. 5. Typical density of states (DOS) calculated with the orbital F and AF ordering at 
x = and 0.3. 

The staggered ordering is therefore stabilized with increasing Jaf (x) in the moderately 
doped region. In the heavily doped region (IV), the expression of Jf,af does not hold any 
more because t ~ [3 e g. There the staggered ordering is unstable due to the lower bandwidth 
than that of the uniform one, leading to the orbital F ordering. 

The competition between the real and the complex orbital states is understood as follows. 
The complex state is stabilized only in the moderately doped region with the advantage of 
the isotropic band structure. In the other region, some other mechanism is rather important 
than the isotropy: In the small doping region, the real state realized in Fig. [I] and |3| 
(x < 0.25) is found to be stabilized mainly due to the hybridization between the occupied 
and the unoccupied bands via the off-diagonal hopping integrals. In the heavily doped region, 
on the other hand, the low dimensional band structure, 9i = On = (two dimensional) or 
7r (quasi-one dimensional), is prefered where the isospin moment is along the z axis (real 
orbital state). This is due to the relative location between the van Hove singularity of DOS 
and the fermi level .0 The fermi level with small electron concentration (x ~ 1, heavily doped 
region) is located near the band edge. The low dimensional band structure with singularities 
at the top and the bottom of the band can therefore lower the kinetic energy effectively with 
large accomodation at the singularity near the fermi level. 

With increasing $/t, the staggered state becomes unstable because Jaf goes to zero 
whereas Jf remains to be finite. This corresponds to the recovery of the orbital Nagaoka 
F. One can therefore expect the uniform complex state in the small doped region with 
the strong correlation limit. The obtained ordering is however the real uniform one with 
\x 2 — y 2 ), not the complex one. The schematic phase diagram in this limit is given in Fig. 
| 
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FIG. 6. Schematic phase diagram in the strong correlation limit. 

This can be understood as follows. In this limit, only the DOS near the band edge matters 
because the bandwidth of the uniform state does not depend on the orbital shape. The DOS 
arises at the edge most sharply with \x 2 — y 2 giving the largest kinetic energy gain. Orbital 
Nagaoka F is therefore realized with \x 2 — y 2 ) in the strong correlation limit. (It should not 
be confused with the real state in the small doping region with realistic parameters (Fig. 
|3|), where the ordering is orbital AF stabilized due to the inter-band hybridization.) 

Curves obtained in Fig. [I] are non-monotonic with a common tangential line contacting at 
two different x, where the curve is convex upwards. This means that the phase coexistence 
with two different concentrations has higher energy than the single phase. Therefore a 
spontaneous phase separation! does not occur in our results. 

The e 9 state specified with the isospin orientation T stabilizes the Jahn- Teller (JT) 
deformation expressed as (T^ ■ Q u + (T^ ■ Q v , where Q u and Q v denote the normal 
coordinates of the displacement of the oxygen ions A a (a = x, y, 

2A 2 - A x - A v A x - A v 

The complex state realized with 0.35 < x < 0.45 corresponds to T//e y . With this ordering, 
therefore, the JT distortion does not occur. The observed disappearance of the JT distortion 
in the spin F metallic regionll!@ might be explained by this type of the orbital ordering. 
Another theoretical proposal is the orbital liquid state where the planer orbitals, \x 2 — y 2 ), 
\y 2 — z 2 ), and \z 2 — x 2 ) are resonating to form a quantum liquid state.! With this resonance, 
the local isotropy is recovered, and the JT distortion disappears on average. The complex 
state can actually be expressed in the form of such a resonance as, 
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_L | x 2 _ ^ ± • 1 | 3z 2 - r 2 ) = ^ [|x 2 - y 2 ) + e ±i2 ^ 3 |^ 2 - x 2 ) + e^/ 3 |y 2 - z 2 ) . 

(8) 

This can be regarded as a formation of the T y -component by the resonance via the transverse 
components = T x + iT y . From Eq. ([8]), the complex orbital ordering corresponds to 
the coherent (in time and space) resonance among the planer orbitals, with the relative 
phases e ±27r//3 being fixed. The orbital liquid state, on the other hand, corresponds to the 
resonance without phase coherence. In this sense, the complex orbital state obtained here 
is the meanfield state akin to the orbital liquid state, and may provide a rough estimation 
of the energy of the latter stateB Because both states give no JT distortion on average, 
direct observations of the orbital state is needed to distinguish them, not via the lattice 
deformation, but via the difference of the spatial orbital correlations. Several probes are 
available, such as the anomalous X-ray scatteringJEH! the X-ray charge density study by 
using of the maximum entropy method (MEM),0 the magnetic Compton scattering,^! and 
the polarized neutron scattering^. 

In summary, we studied the competitions among the uniform complex, staggered com- 
plex, and real orbital states in CMR compounds. In the moderately doped region, the com- 
plex orbital state is stabilized, where the ordering changes from the canted one (x = 0.25, 
0.3) to the staggered one (0.35 < x < 0.45). This can be understood in terms of the com- 
petition among the band narrowing and the gap associating with the staggered structure, 
and the hybridization with the unoccupied band. The staggered complex ordering is not 
accompanied with the Jahn- Teller deformation. The obtained complex orbital is a coherent 
resonance among the planer orbitals with constant phase angles. If the coherency is lost, 
the state reduces to the orbital liquid state which can be distinguished by the observation 
of the spatial orbital correlations. The phase separation does not occur with the complex 
orbital state obtained here. 
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